Low-temperature magneto-thermal transport investigation of a Ni-based 
superconductor BaNi 2 As 2 : Evidence for fully gapped superconductivity 
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We have performed low-temperature specific heat and thermal conductivity measurements of the 
Ni-based superconductor BaNi2 As2 (T c = 0.7 K) in magnetic field. In zero field, thermal conductivity 
shows T-linear behavior in the normal state and exhibits a BCS-like exponential decrease below 
T c . The field dependence of the residual thermal conductivity extrapolated to zero temperature is 
indicative of a fully gapped superconductor. This conclusion is supported by the analysis of the 
specific heat data, which are well fit by the BCS temperature dependence from T c down to the 
lowest temperature of 0.1 K. 

PACS numbers: 74.70.Dd,74.25.Fy,74.25.Op 



Since the discovery of superconductivity in 
LaFeAs(0,F) Q, there has been considerable inter- 
est in the oxy-pnictide RTPn(0,F) and the related 
structure AT 2 Pn 2 (i? = La, Ce, Sm, Nd, yl = Ca, Ba, 
Sr, Eu, T = Fe, Ni, Pn = P, As) due to (i) their high 
superconducting temperature T c (up to 55 K for Sm- 
FeAs(O,F)01 and 38 K for (Ba,K)Fe 2 As 2 Q), (ii) their 
proximity to magnetism, and (iii) the large variety in 
structure and composition that supports superconduc- 
tivity. Ni-based materials differ from their Fe-based 
cousins in that (i) long-range magnetic order has not yet 
been observed in close proximity to superconductivity, 
although similar structural transitions are found, and 
(ii) T c does not exceed 5K in any of the Ni-based 
systems, although for virtually every structural variant 
where the Fe system superconducts, so too does the Ni 
analog. Substitution of Ni for Fe modifies a number of 
properties that may be responsible for suppression of 
T c . Identifying the commonalities and, especially, the 
differences between the compound in the two families, 
will in turn provide clues into why T c in Fe-based 
superconductors is so high. 

One important issue to resolve with any new supercon- 
ductor is to identify the superconducting order parame- 
ter as that may shed light into the pairing mechanism. 
Comparing the character of the superconducting gap in 
Ni- and Fe-based systems present an ideal opportunity to 
elucidate the pairing mechanism responsible for the rela- 
tively high T c in Fe-based compounds. In the iron-based 
oxy-pnictides or AFC2P712 compounds, NMR measure- 
ments suggest a nodal gap structure 0, Q , which con- 
tradicts penetration depth [|| , Andreev spectroscopy Q , 
and angle-resolved photoelectron spectroscopy (ARPES) 
measurements Q that indicate a fully gapped multi- 
band superconductor. To our knowledge, similar mea- 
surements that address the character of the supercon- 
ducting order parameter have not been done on the Ni 
analogs [|, M, [n|, (H 0, 0, [H] . 



We present a study of the gap structure of BaNi2As2 
using low temperature thermal conductivity and heat ca- 
pacity. Thermal conductivity, as a bulk probe of low 
energy delocalized excitations, has proven to be a power- 
ful tool for identifying nodal [3l , s-wave [l?} , and multi- 
band superconductors [3]. We chose to study BaNi 2 As2 
as a representative Ni-based superconductor. Bulk su- 
perconductivity in this material at 0.7 K is confirmed 
by heat capacity measurements, while, similar to many 
Fe-pnictides, a first order structural transition exists at 
13OK0. 

Single crystals of BaNi 2 As2 (ThCr 2 Si2-type tetrag- 
onal structure) were grown in Pb flux as described 
in Ref. 12. Thermal conductivity was measured be- 
tween 40 mK and 4K by a standard one-heater and 
two-thermometers technique on plate-like crystals with 
dimensions of ~lx0.5x0.lmm 3 for a heat current 
q || [100]. Specific heat of several single crystals (6.91 mg) 
were measured by a standard heat pulse method. In the 
thermal conductivity study, two samples from different 
batches were measured, which are denoted as #1 and 
#2. Sample #1 is the same one used in previous resis- 
tivity experiments fRcf.[l2l). Magnetic field was applied 
within afo-plane for HJ^q (#1 and #2) and H\\q (#1). 

First, we consider the possible superconducting gap 
symmetry of BaNi 2 As 2 based on heat capacity data. Fig- 
ure QJ a) shows specific heat divided by temperature C/T 
as a function of T in zero field which supports the notion 
that BaNi 2 As 2 is a weak coupling superconductor. The 
upturn below 0.2 K is ascribed to a nuclear quadrupolar 
Schottky anomaly arising mainly from As. The solid line 
is a fit to C = Cbcs + C Sc h, where C Sch =A/T 2 and C* B cs 
is the BCS expression: 
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Here t = T/T c , / = l/[cxp(£/fc B T)+lh £ = (e^A 2 ) 1 / 2 
and y — e/Ao, as described in Rcf.lig. As seen in 
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FIG. 1: (Color online) (a) Zero-field low-temperature specific 
heat, plotted as C/T vs T, of BaNi 2 As 2 and Sr 2 Ru0 4 [H[ 
(with T c and *yT c scaled to that of BaNi 2 As 2 ) for comparison. 
The solid curve is a fit to the theoretical expectation for BCS 
superconductivity including a nuclear quadrupolar Schottky 
contribution, (b) C/T vs T in fields of 0, 0.01, 0.02, 0.03, 
0.04, 0.05 and 0.1 T (right to left as indicated by the arrow) 
for H\\ab. (c) Upper critical field H C 2 for H\\ab as determined 
by the midpoint of the jump in C/T. The solid line represents 
a least-squares fit to the data. 



Fig. [TJa), the specific heat data can be fit over a 
wide temperature range from just above T c down to 
the lowest temperature, yielding a Sommerfeld co- 
efficient 7 = 12.3mJ/molK 2 , A = 8.6 x 10" 3 mJ K/mol, 
equivalent to the energy splitting of the nuclear 
quadrupolar levels of As of i/q « 30 MHz [gcj , and 
A = 0.0946 meV = 1.61 k B T c (T c = 0.68 K was fixed dur- 
ing the fit). The reduced magnitude of the en- 
ergy gap Ao, compared with the expected BCS gap 
Abcs =0.103meV= 1.76 fee T for a weak coupling BCS 
superconductor, is likely due to the fact that our 
BaNi2As2 samples are in the dirty limit as will be dis- 
cussed below based on the analysis of thermal conductiv- 
ity data. In contrast, specific heat of a well established 
unconventional superconductor Sr2Ru04, also displayed 
in Fig. Ufa), shows a much slower decrease with temper- 
ature, expected for nodal superconductivity [2l[ . 

Fig. QJb) shows the temperature dependence of C/T 
in several fields for H\\ab in BaNi2As2. As field in- 
creases, the jump shifts to low temperature region, 
whereas it becomes broader and less detectable above 
0.05 T. The upper critical field if C 2, determined by 
the midpoint of the jump at each field, is displayed 
in Fig. []Jc). The solid line is a least-squares fit 
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FIG. 2: (Color online) Temperature dependence of the ther- 
mal conductivity k(T) in BaNi 2 As 2 for heat current q||[100], 
in zero field and in the normal state (0.5 T) for H\\ab. Arrows 
indicate T c determined from the specific heat measurement. 
Solid square and cross symbol represent the electronic thermal 
conductivity K e = L T/p with L = 2.45 x 10~ 8 Wfi/K 2 for 
and 0.5 T, respectively, derived from resistivity data using the 
Wiedemann-Franz law. Dotted line in main figure and solid 
line in the inset show a fit to k(T) with a BCS curve defined 
as K = Cexp(-aT c /T). 



to the H C 2 data, yielding an initial slope dH C 2/dT c 
=— 0.226 T/K, which results in a zero-temperature up- 
per critical field H C 2(0) = 0.11 T, using the estimate 
H c2 (0) = -0.7T c dH c2 /dT c This yields a Ginzburg- 
Landau coherence length £ = 550 A from the relationship 
t=($o/2irH c2 (0)) 1/2 , where <J> = 2.07 x 10~ 7 Oecm 2 is 
the flux quantum. An estimate of the electronic mean 
free path l e is obtained from the normal state thermal 
conductivity using K/T—l/3jv-pl c , with a Fermi veloc- 
ity up = 2.96 x 10 5 m/s which was calculated by methods 
described in RefS We find Z C «70A and f e /f = 0.13, 
which places BaNi2As2 in the dirty limit. The field de- 
pendence of the residual linear term 70T of the heat ca- 
pacity can be used to identify whether or not the super- 
conducting order parameter has nodes. However, in the 
low temperature region, especially below 0.2 K, the ad- 
denda contribution, as well as the nuclear term, grow 
with increasing field, making an accurate estimate of 
7o(-ff) difficult. In this regard, thermal conductivity, 
which is not subject to these effects, is ideally suited for 
directly determining the low-temperature behavior. 

Fig. [2] shows the temperature dependence of thermal 
conductivity n(T) of BaNi2As2 in zero field and 0.5 T 
(^>H C 2) for heat current g||[100]. In zero field and in the 
normal state above T c = 0.7K, k(T) exhibits an approxi- 
mately T-linear variation, which continues to the lowest 
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temperature in a field of 0.5 T. Using the Wiedemann- 
Franz law, we estimate an electronic thermal conductiv- 
ity in the normal state K e = L T/p in zero-field and 0.5 T 
where the Lorenz number L = 2.44x 10~ 8 Wf2/K 2 . As 
seen in Fig. [21 K e (T) is very close to k(T) in the nor- 
mal state, suggesting that heat transport in the nor- 
mal state is dominated by the electrons. In the su- 
perconducting state, k(T) follows an exponential form 
(~ exp(— aT c /T) with a = 1.34) down to 0.2 K as shown 
in the inset of Fig.[2l This evolution of k(T) in BaNiaAsa 
is quite similar to the conventional superconductors tin, 
aluminum, and zinc with T c = 3.7, 1.2 and 0.84K, re- 
spectively (24| ■ The thermal conductivity for these three 
elements is also described by k(T) ~ exp(— aT c /T) with 
a = 1.3- 1.5 below T c . In fact, the thermal conduc- 
tivity based on BCS superconductivity for weak cou- 
pling [25| is well-described by the formula exp(— aT c /T) 
with a =1.47. Thus, the a-value of 1.34 for BaNi2As2 
is in good agreement with that for conventional BCS 
superconductors. Below 0.2 K, «(T) does not fall as 
rapidly as expected for an exponential dependence and 
reaches 2xlO _4 W/Km at 50 mK. Similar deviations 
from an exponential dependence were observed in a num- 
ber of conventional superconductors, including tin and 
aluminum [3] , and are usually attributed to phonon con- 
tribution. When most of the normal quasiparticles are 
frozen out, the intrinsic electronic thermal conductivity 
is too low to make an appreciable contribution. In fact, 
an upper limit of the phonon thermal conductivity esti- 
mated via Kp h = i/3T 3 {u)/ ph = 2.0xlO- 4 W/Km, using 
the phonon specific coefficient (3 = 18.7 J/K 4 m 3 , mean 
phonon velocity (v) = 1860 m/s |1 2T and the phonon mean 
free path Z p h = 1.38 x 10 -4 m [26J, is equal to the exper- 
imental value 2x 10 _4 W/Km at 50 mK in BaNi 2 As 2 . 
Although we cannot rule out other origins of the low 
temperature tail, agreement between the estimated and 
the measured values at 50 mK suggests a phonon origin of 
the low temperature thermal conductivity in BaNi2As2- 
Note that nuclear quadrupolar moments, while leading 
to a Schottky term dominating specific heat below 200 
mK, do not contribute to thermal conductivity. 

The low-temperature magnetic field dependence of k 
is instructive for determining the gap structure. Fig. [3] 
shows the low-temperature expansion of k/T vs T 2 of 
BaNi2As2 for a heat current g||[100] in several fields for 
H±q. The straight lines are fits to k/T = k /T + bT 2 , 
where kq /T is the residual term extrapolated to T = K 
at each field. With increasing field, kq/T rapidly in- 
creases above 0.02 T and saturates above 0.15 T. This 
rise could only be attributed to the electron contribu- 
tion, as the phonon thermal conductivity may only go 
down in magnetic field due to additional scattering from 
vortices in the mixed state. The field dependence of the 
residual thermal conductivity of BaNi2As2 is shown in 
Fig. Ha) as (k /T)/(k„/T) vs H/H c2 for H±q {*!). 
From the thermal conductivity data, H C 2 = 0.16T is dc- 



0.20 



0.15 - 



I 

H 

2 



0.10 



0.05 - 



0.00 



1 r 

BaNi 2 As 2 (#1) 
100], H±q 




0.20T 
0.15T 
0.12T 
0.10T 
0.09T 
0.08T 
0.07T 
0.06T 
0.05T 
0.04T 
0.03T 
0.02T 
0.01T 
0T 



0.00 



0.01 



0.02 

2 2 

T(K) 



0.03 



0.04 



FIG. 3: (Color online) Low-temperature dependence of k/T 
vs T 2 for q\\ [100] in several fields for H.Lq. The solid lines are 
fits to the data in the temperature range by k/T = kq /T + bT 2 
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FIG. 4: (Color online) (a) Residual thermal conductivity 
ko/T of BaNi2As2, normalized by the normal state value 
K n /T above H C 2, as a function of H/H C 2 for H.Lq ( # 1). 
A small contribution from the Pb flux has been subtracted. 
For comparison, data for several superconductors with dif- 
ferent superconducting gap characteristics are shown: Nb 
(clean, fully gapped s-wave) J22|, InBi (dirty, fully gapped s- 
wave) [13], MgB 2 (multi-band gap) 0, Tl-2201 (d-wave with 
line nodes) [la]. The dotted lines are guides to the eyes, (b) 
{ko/T)/(k u /T) vs H/H c2 of BaNi 2 As 2 for H±q (*1*2) and 
H\\q {*!). 
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tcrmined as the field where kq/T becomes independent 
of the magnetic field. n n /T at 0.2 T is used as the nor- 
mal state value. Data for several superconductors with 
different gap structures are also shown for comparison. 
The simplest test of s-wave superconductivity is whether 
the residual thermal conductivity exhibits a concave vari- 
ation with field as vortices begin to penetrate the sam- 
ple at low fields. The exponential- like concave behav- 
ior is understood as follows: when vortices first enter 
the sample at H>H C \, quasiparticles that contribute to 
thermal conductivity are mostly localized around the vor- 
tex cores. Therefore, at low field, those quasiparticles 
cannot carry heat until the intervortex spacing is de- 
creased sufficiently that states created in the vortices be- 
gin to overlap, become delocalized, and hence, transport 
heat, as illustrated by the clean s-wave superconductor 
Nb[27|). In contrast, the presence of nodal quasiparticles 
leads to a convex field dependence as a consequence of 
the Volovik effect [H], as illustrated by Tl-2201,0. As 
seen in Fig. H^b), (kq/T)/(k d /T) exhibits concave fea- 
ture independent of heat current direction with respect 
to magnetic field, and independent of sample difference. 
Consequently, the concave field- dependence of the ther- 
mal conductivity implies that BaNiiAs2 is a fully gapped 
superconductor. 

The field dependence of n(H) for BaNi2As2 more 
closely resembles the dirty s-wave superconductors such 
as InBi, than the clean case (i.e., Nb), which is consis- 
tent with our estimate of l e /^ <1. In addition, the data 
appear to display a shoulder-like anomaly at 0.7 H C 2, sug- 
gesting an additional energy scale at that field. This may 
reflect a spread in H C 2 as reflected by the width of the 
heat capacity anomaly in field. On the other hand, the 
shoulder might be due to multiband character of super- 
conductivity in BaNi2As2, although the anomaly in ther- 
mal conductivity is much less dramatic than in MgB2. 
And there is no obvious evidence for multiband physics 
in the heat capacity data of Fig. [1] However, interband 
scattering would be expected to wipe out multiband ef- 
fects in our relatively dirty material. Additional studies 
on cleaner crystals would help to resolve this issue. 

For the FeAs-based superconductors, it has been ar- 
gued that an unconventional superconducting state may 
still give rise to a fully gapped excitation spectrum if the 
nodal planes, where the phase of the order par ameter 
changes sign, do not intersect the Fermi surface [29j . Den- 
sity functional calculations for BaNi2As2, however, find 
a much more complicated Fermi surface for which it is 
difficult to imagine finding a nodal plane which does not 
intersect the Fermi surface 23(. Thus, the fully gapped 
Fermi surface in BaNi 2 As2 is due to the superconduct- 
ing order parameter that does not change sign. If the 
superconducting order parameter in FeAs systems does 
indeed change sign, the pairing mechanism in the FeAs 
and NiAs systems is likely to be of different origin, as 
argued by band structure calculations 23|, 3p| ] . 



In conclusion, we have performed magneto-thermal 
transport experiments to elucidate the superconducting 
gap symmetry of BaNi2As2. The following results indi- 
cate that BaNi2As2 is a dirty fully gapped superconduc- 
tor: (i) specific heat data are well-fit by the BCS for- 
mula, (ii) the coherence length is much larger than the 
electronic mean free path, and, (iii) the field dependence 
of the residual thermal conductivity is consistent with a 
dirty fully gapped superconductivity. 
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